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C> 1 . Introduction which are observed in the optical in diffuse clouds. 

i~ I ' The Herschel Space Telescope permits observations of hy- 

; Planets are formed in discs around young stars. The building drides and their ions (OH , H+, CH, CH + , ...), and warm 
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ABSTRACT 

Despite its importance in the thermal balance of the gas and in the determination of primeval planetary atmospheres, the chemistry 
in protoplanetary discs remains poorly constrained with only a handful of detected species. We observed the emission from the disc 
around the Herbig Be star HD 100546 with the PACS instrument in the spectroscopic mode on board the Herschel Space Telescope 
as part of the GaS in Protoplanetary Systems (GASPS) programme and used archival data from the DIGIT programme to search 
for the rotational emission of CH + . We detected in both datasets an emission line centred at 72.16 /jm that most likely corresponds 
to the .7=5-4 rotational emission of CH + . The 7=3-2 and 6-5 transitions are also detected albeit with lower confidence. Other CH + 
rotational lines in the PACS observations are blended with water lines. A rotational diagram analysis shows that the CH + gas is warm 
at 323+ 2 j 20 K with a mass of ~3 x 10~ 14 -5 x 10~ 12 M Q . We modelled the CH + chemistry with the chemo-physical code ProDiMo 
using a disc density structure and grain parameters that match continuum observations and near- and mid-infrared interferometric 
data. The model suggests that CH + is most abundant at the location of the disc rim at 10-13 AU from the star where the gas is warm, 
which is consistent with previous observations of hot CO gas emission. 
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of planets is intimately linke d to the grain-grow th process and light mo i eC ular gas (e.g., H 2 0), which all emit in the far- 
the disc's chemical evolution (|Bergin et alj2007|). The chemical infrared. Rotational lines of CH + at 300-500 K in the far- 



and temperature profiles of discs determine the chemical com- IR have been observe d by ISO towards the NGC7027 PDR 

H ' position of the gas that will eventually be incorporated into a dCernicharo et al.ll997l) . The isotopologue 13 CH + 7=1-0 was de- 

9 .' planet's atmosphere, once a solid ~ 10 Earth mass core is formed tected in absorption in the diffuse medium by iFalgarone et ali 

dArmitage||2010|). Millim etre observations of molecules in discs (Fjoosh . The lowest rot ational transition was detected in emis- 

probe the outer cold p art (iDutrev et al.ll997tlvan Zadelhoff et all sion in the Qrion Bar dNavlor et al.l l2010[ lHabart et al.l|20Tob 



2001; Thi et al. 2004j|Semenov et alj2005|;|Henning et al.|2010; and in absorption towards DR 21 (Fakarone et al. 2010) with 
Oberg etalj |2010l e.g.). Minor species located in the inner Herschel. The neutral counterpart CH has also been detected 



disc, such as C ? H? and HCN, were observed with the Spitzer in me diffuse medium by Herschel (lOin etal.ll2010l iGerin et al.l 
Space Telescope (Carr & Naiita 2008; Pascucci & Sterzik||2009[ Fj()ToT) . In this letter, we discuss the detection of emission lines 
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lLahuis et alj|2006|; |Pontoppidan et al.||2010|). Gas-phase carbon with Herschel that we assign to the rotational transition of 

is believed to be locked into C + , C, and CO. However, a sig- CH + from the disc around the 10 Myrs old Herbig Be star 

nificant fraction of the carbon can be locked into CH and CH + , H D 100 546 (M*= 2.2 M , L»= 27 L , r eff = 10500 K) located at 

103 pc (Ivan den Ancker et al.lFl997h . HD 100546 shows a very 

* Herschel is an ESA space observatory with science instruments pro- rich gas and a solid spectrum in the infrared and millimetre 
vided by Principal Investigator consortia. It is open for proposals for 
observing time from the worldwide astronomical community. 
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Fig. 1. Continuum-subtracted Herschel-PACS spectra of HD 100546 around the CH + 7=6-5, 4-3, 3-2, 2-1 (DIGIT data), and 7=5- 
4 (GASPS and DIGIT data) lines. The 3<x statistical error levels do not include the 30% calibration uncertainty. A rotational diagram 
using those transitions is plotted in the lower-right panel. The errors in the diagram are the quadratic sum of all the statistical and 
calibration errors (err tot = ~<J(3cr stat ) 2 + (0.3F„i, s ) 2 ). Blended lines are considered as upper limits (red filled dots). 



range s (Brittai n et al]|2009t ISturm et al.ll2010t|Panic et al.ll2010t 
Malf ait etalJll998l) . 



2. Herschel observations and rotational diagram 



We observed HD 100546 with PACS (Poglitsch et al. 2010) on 
board the Herschel Space Telescope dPilbratt et al.ll2010h in the 
72-73 yum spectroscopic scan range (Av= 163 km s~'), as part of 
the open-time key programme GASPS (GaS in Protoplanetary 
Systems, obsid 1342188437), which aims at a systematic study 
of gas and dust in discs around stars aged between 0.5 and 30 
Myrs. The observations were reduced using the standard pipeline 
in HIPE 4.2. We detected a line at 72.14 yum that we assign 
to the rotational 7=5-4 emission of CH + (see Fig. [TJ. We also 
detected the same line in the s cience demonstra tion data spec- 
trum taken by the DIGIT team dSturm et al.f2 010). together with 
the 7=6-5 (60.24 yum) and 7=3-2 (119.68 fim) lines but with a 
lower signal-to-noise ratio. It was not possible to strengthen our 
analysis with detection of either the lowest rotational transition 
7=1-0, which lies in the SPIRE/HIFI range, or the 7=2-1 and 
7=4-3 lines in the DIGIT dataset, which are blended with water 
emission lines at the PACS resolution. The ortho H2O 2i2-loi 



emission at 179.52 fim is blended with the CH + 7=2-1 line at 
179.61 yum, and the para-H20 322-2n at 89.98 yum is blended 
with the CH + 7=4-3 line at 90.02 //m. The flux calibration accu- 
racy of the Herschel-PACS instrument is currently 30% in ad- 
dition to the statistical noise. In addition, flux uncertainties are 
introduced by difficulties in determining the continuum level in 
the DIGIT spectra. A new reduction with HIPE 7.0.822, which 
includes division by the spectral response function and flatfield- 
ing, does not improve the noise level in the continuum. The 
line ratio [O I] 63 /mi/CH + 7=5-4 is 80. We fitted a Gaussian 
profil e to each line. The observed CH + [O I], [C II] dSturm et all 
1201 Ol) . and CO 3-2 dPanic et alj|2010b fluxes are all reported in 
TableQ] 

We constructed a rotational diagram by assuming LTE 
population and optically thin emission (see lower right panel 
of Fig. [T). The level energies, line frequencies, and Einstein 
sponta neous emission coefficients were computed by Miiller 
(2010) based on the spectroscopic analysis of Amano (2010) 
. We analysed the diagram using censored data regression 
and bootstrapping methods. We derived an excitation temper- 
ature of 323^j5Q° K, a disc-averaged CH + column density of 
4.3^3 3 j 8 x 10 12 cirT 2 , consistent with the upper limit found in 
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Table 1. Observed and modelled line fluxes 



Transition 


A 


obs. 


model 1 


model 2 




(/urn) 


(10" 


17 \\r -2\ 

W m ) 




[0 I] 3 P!- 3 P 2 


63.19 


554.37 ± 167° 


785 


763 


[O I] 3 P - 3 Pi 


145.54 


35.70 ±11.4 


35.7 


42.7 


[C II] 2 P 3/2 - 2 P 1/2 


157.75 


31.87 ± 10.0 


12.7 


17.2 


p-H 2 3 22 - 2n 


89.90 


< 14.32 ± 5.6* 


4.6 


7.4 


CH + J = 4 - 3 


90.02 


< 14.32 ± 5.6* 


3.7 


6.3 


CH+ J = 6 - 5 


60.25 


10.32 ± 5.7 


4.2 


8.9 


CH + J = 5 - 4 C 


72.14 


6.86 ± 3.4 


3.3 


6.7 


CH+ 7 = 3-2 


119.87 


2.16 ± 1.1 


1.8 


2.9 


CO 7 = 3 - 2 d 


866.96 


0.10 ±0.03 


0.08 


0.24 



8 The total (3cr statistical + 30% calibration) errors are given for the 
PACS observations. * The blended H 2 0+CH + lin e is detected c The 
PACS values are from the DIGIT programme (ISturm et al.l l2010h 
except for the CH + 7 = 5-4 flux (GASPS programme). d Data from 
iPanic etal.U20ld) with 3<x error. 



Table 2. Model parameters. 





MCFOST" 








inner disc 


surf, layer 


outer disc 


Inner radius R m (AU) 


0.24 


13 


13 


Outer radius R out (AU) 


4 


50 


500 


Surf, density exponent q 


1 


0.5 


1.125 


Scale height TZiooau (AU) 


6 


14" 


14° 


Scale height exponent /3 


1.0 


1.0 


1.0 


Total dust mass M dusC (M Q ) 


1.75(-10)* 


3(-7) 


4.3(-4) 


Dust mass (a < 1mm, M Q ) 


1.75(-10) 


3(-7) 


1.3C-4) 


Min. grain radius a m j„ (pm) 


0.1 


0.05 


1 


Max. grain radius a max (jum) 


5 


1 


10 4 


Grain power law index p 


3.5 


3.5 


3.5 


Silicate grain density (g cirr 


3 ) 3.0 


3.0 


3.0 




ProDiMo^ 






ISM UV field (x, Habing) 




1.0 




viscosity (a) 




0.0 




Non-thermal speed (km s~') 




0.15 




Disc inclination (°) 




42 




UV excess 




0.013 




UV power-law index 




6.5 




Cosmic ray flux <r(s~') 




1(-17) 




PAH C 15 oH 3 o mass (M Q ) 




1.8C-7) 




Gas mass (M©) 




5(-4) (model 1) 




Gas mass (M Q ) 




l(-3) (model 2) 





three parts: the inner disc, the outer disc, and an upper layer on 
the top of the outer disc (see Table[2] Fig. IA. II and the sketch in 
iBenisty et al.ll2010|). The inner and outer discs a re separated by 
a gap dBouwmanetalJl2003l:lGradvetal.ll2005h . The fit to the 
SED in Tatulli et al. (2011) suggests a scale height value of 10 
AU at 100 AU, which is consistent with the scattered-light im- 
ages, but equally good fits can also be obtained with other values. 
Therefore we tried several disc structures by varying the scale 
height Hiooau- Using the disc structures derived from the fit to 
the SED, we modelled the gas chemistry and radiative transfer in 
the HD 100546 di sc with the ProDiMo code dWoifke et alj2 009: 
Kamp et al] |2010h . The ProDiMo code computes the chemistry 
and heating-cooling balance of the gas self-consistently. We as- 
sumed a non-viscous gas (ff=0) and a PAH mass (C1 50H30) that is 
consi stent with the strong PAH emissions in the IR dKeller et al.l 
2008). Large PAHs (> 100 car bon atoms) exhibit a stro ng fea- 
ture between 10.9 and 11.3 fim dBauschlicher et al.l2009l) as seen 
in the spectrum of HD 100546. The amount of PAH influences 
both the thermal state of the disc and the ionisation state of the 
upper layers. The non-thermal line width is constrained by the 
resolved CO 3-2 profile. We used a deredde ned FUSE and IUE 
UV s pectrum as input to the gas modelling (Martin-Zaidi et all 
2008). We also ran a ProDiMo model where the disc vertical 
structure is calculated according to the gas pressure and found 
that the height of the rim at 13 AU is consistent with the fixed 
three-part structure. 

The main CH + formation reaction C + + H 2 -> CH + + H 
(reaction 1) is endothermic by 0.398 eV (4537K). The re- 
action is efficient only at a few 100 K found in shocked 
or turbulent regi o ns or when H2 is vibrationally excited 
dHierl et all Il997t lAgundez et all l2010h . The rate for reac- 



tion 1 in iHerbst & Knudsonl dl98lh underestimates the exper 



8 Values taken from lBenistv et all ( 120101) and Tatulli et al. (201 1) except 
for the scale height H loaAV at 100 AU. * a(-fi) means a x 10^. c This 
work. 



imental data by a f actor of 2-3, and the UMIST 2006 rate 
(Woodall et a n i2007l) is a factor 7-8 lower than the experimen- 
tal values. Initial models with a low rate for reaction 1 un- 
derpredict the CH + abundance b y more than a fac tor 10. We 
adopted the value suggested bv iHierl et al.l (1199 7): k\(T) = 
(7.4 + 0.8) x 10" 10 e" 4537/r cm 3 s _1 . We also take the enhanced 
rate for vibrationally-excited molecular hydrogen into account 
by assuming that the r ate is equal to the Lang evin collision value 
of 1.6 xlO" 9 cm 3 s" 1 dAgtlndez et alj |2010). The alternative ra- 
diative association reaction C + +H — » C H + +hv has an extremely 
low rate of ~ 10 16 -10~ 17 cm 3 s 1 (Ba rinovs & van Hemertl 
2006). The warm and UV-excited gas at the disc atmosphere 
around a Herbig AeBe star renders reaction 1 fast. CH + reacts 
with H2 to form CH 4 ,, which in turn reacts with H2 to form 
CHo . As both exothermic reactions are fast, an enhanced abun- 
dance of CH + is accompanied by a high abundance of CHJ and 
observing CH + at 4232.7 A (Martin-Zaidi et al., submitted), and CH+, which have not been detected yet. CH + can also dissocia- 
a lower limit to the CH + mass of 3.3+f 9 x 10~ 13 M Q . The CH + tively recombine with an electron into C + and H or exchange 

its charge with PAHs. CH + is also destroyed upon absorption 
of UV photons. The absorption cross-section peaks at ~950 A. 
Therefore, the photodissociation around a B9.5V star without 
strong UV excess due to accretion is not as fast as with an inter- 
stellar UV field. The photodissociation and photoionisation rates 
are computed inside the disc using cross-sections taken from 
van Dishoeck et al. (2006). The steady-state chemistry includes 
188 gas and solid species. 

We ran several models with disc gas mass between 10~ 4 
and 10~ 2 M G and Hiooau between 10 and 20 AU (Fig. lATll 
In Table [TJ we show the three detected line fluxes for CH + out 
of the 15 rotational transitions in the model, the OI and CII 
fine-structure lines, and the CO 3-2 transition for our preferred 
models, which have a disc mass of (0.5-1) x 10~ 3 M and 



emitting gas is warm, which suggests that CH + is located within 
100 AU of the star. The CH + abundance is > 3.9 x 10" 9 rela- 
tive to an upper limit of warm molecular gas from the VISIR 
H? S(l) 17 urn observati on of M(7g as =300 K)< 8.4 x 10" 5 M Q 
dMartin-Zaidi et al.ll2010l) . 

3. Modelling the physics and chemistry of CH 

We first fitted the spectral energy distribution (SED) and infrared 
interferometric VLT-AMBER and MIDI data to constrain the 
dust properties and the disc structure with a hydrostatic disc 
using the MCFOST radiative transfer code (Pinte et al. 2006, 
2009). The details of the fit are discussed in Benistv et al.ld2010l) 
and Tatulli et al. (2011, submitted). The disc is composed of 
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1 10 100 



Fig. 2. CH + abundance structure in the HD 100546 disc with 
respect to the number of H-nuclei. The H2 abundance relative to 
total H-nuclei abundance contours in log-scale are overplotted 
in blue, while the r gas =300 and 3000K contours are overplotted 
in white. The line-of-sight CH + observations in the UV would 
cross the disc at z/i-0.9 (42 °). 



Hiooau = 14 AU, for the outer disc. All the modelled lines agree 
within a factor 2 with the observations. In particular, the sum of 
the P-H2O 322 - 2n and CH + 7 = 4-3 fluxes is consistent with 
the observed valu e. The col l ision rates of CH + with electrons 
were computed bv lLim et al.l (1 19991) . In the absence of published 
collision rates with H2, we scaled the collision rates with He by 
1.39 for J < 6 transitions (iHammami et al J 120091: iTurpin et al.l 
12011 . 



4. Discussion 

The models suggest that CH + is mostly located on the rim (10- 
13 AU) of the outer disc and at the disc surface (Pig . |2j . Although 
CH + is abundant in the tenuous inner disc, its amount is insuf- 
ficient to contribute to the 72.14 fim emission. CH + is located 
where the gas is mostly heated by collisional de-excitations of 
nascent and UV-pumped HJ. Part of the CH + is formed by the 
reaction with excited H* v but the main fraction comes from the 
reaction of C + with thermal H 2 at T m > 400 K. CH + reaches 
a maximum abundance of ~ 10~ 7 just behind the rim and has a 
total mass of 10 13 -10 12 M . The disc average CH + abundance 
is (0.7-7) x 10 9 relative to H-nuclei, consistent with the simple 
analysis in Sect. [2] The high CH + abundance in the rim region 
and upper disc layers where hydrogen is in atomic form does not 
affect the C + abundance much, which reaches 10~ 4 . However 
below the disc atmosphere, H2 self-shielding is efficient, and C + 
is converted rapidly into CH + . The formation reaction of CH + 
(C + +H2) competes with the recombination reaction of C + to 

form neutral ca rbon. 

Warm OH (ISturm et al.l|2010h and fluorescent ro- vibr ational 
CO emission dBrittain et alJl2009t Ivan der Plas et al. 2009) have 
also been detected in HD 100546 . OH is mostly formed by the 
reaction of atomic oxygen O with thermally hot H2 and vibra- 



tionally excited H^. The presence of hot CO, CH + and OH sup- 
ports the idea that hot and excited gas chemical reactions occur 
on disc surfaces and at the inner rim. 

5. Conclusions 

We detected the 7=5-4 rotational transition of CH + at 72.16 //m 
in two Herschel datasets. We also tentatively detected other 
lines from CH + at 60.25 and 119.87 fim. Other CH + lines are 
blended with water lines. We modelled the CH + line fluxes us- 
ing the most recent chemical and collisional rates. Searches for 
CH + in other Herbig AeBe discs are warranted to test whether 
the presence of CH + is unique to HD 100546 or if it is related to 
the presence of a high temperature rim of gas, or both. 
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Fig.A.l. The input disc density structure for model 2 con- 
strained by the fit to the SED (see Tatulli et al. 201 1 for details). 
The gap is filled with a gas at density of «h=100 cirT 3 . 
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Appendix A: Disc density and gas temperature 
structure 

We provide for completeness the disc density and the dust and 
gas temperature structures for model 2 (Fig. IA.ll IA.2J and IA.3I ). 
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Fig.A.2. The disc dust temperature profile computed by 
MCFOST for model 2. 
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Fig.A.3. The disc gas temperature profile computed by 
ProDiMo for model 2. 



